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Disorder-induced neutral solitons in degenerate ground state polymers 
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Abstract 

We study the effects of weak off-diagonal disorder on 7r-conjugated polymers with a doubly degenerate ground-state. 
We find that disorder induces a finite density of neutral solitons in the lattice dimerization of a polymer chain. Interchain 
interactions result in a linear potential between the solitons and, if sufficiently strong, bind them into pairs resulting in 
an exponential suppression of the soliton density. As neutral solitons carry spin 1/2, they contribute to the polymer's 
magnetic properties. We calculate the magnetic susceptibility and suggest measurements of the magnetic susceptibility in 
trans-polyacetylene at low temperatures. 

Keywords: Polyacetylene and derivatives, Magnetic measurements, Ising models, Order-disorder phase transitions 



1. Introduction 

Conjugated polymers belong to a large class of quasi- 
one- dimensional Peierls materials, in which the lattice dis- 
torts due to its interaction with itinerant electrons. Much 
attention has been devoted to trans-polyacetylene, which 
has a doubly degenerate ground state. The degeneracy al- 
lows for topological excitations, solitons, which are kinks in 
the lattice dimerization accompanied by a local distortion of 
the electron density. Solitons survive the presence of on-site 
Coulomb repulsion. Although solitons are usually associ- 
ated with excited states, we show in this contribution that 
off-diagonal disorder may induce a finite density of kinks 
in the ground state of trcms-polyacetylene. We discuss the 
effects on the magnetic response. 

2. Mapping on the Random-Field Ising Model 

We describe the polymer chains by a Peierls-Hubbard 
model, which accounts for both a static electron-lattice in- 
teraction and electron-electron interactions. The lattice 
dimerization at position n of the chain is given by: 



A(n) = Ai at {n)+r){n). 



(1) 



Here, Aj a t(n) describes the Peierls distortion and corre- 
sponds to the alternating part of the electron hopping am- 
plitude along the chain direction. The second term, rj(n), 
represents the disorder and stems from the small fluctua- 
tions in this amplitude. While in the absence of disorder 
the ground state energy does not depend on the sign of the 



dimerization A(n) = ±Ao, the weak-disorder correction to 
the energy, 



5E oc - 2J Aj at (n)?j(n) , 



(2) 



removes this degeneracy. As we explained in Ref. [1], the 
energy correction Eq. (El) stabilizes neutral solitons in the 
lattice dimerization of a chain's minimal-energy lattice con- 
figuration. 

As a consequence of the form Eq. (^), the statistics 
of neutral solitons in a weakly disordered Peierls-Hubbard 
chain can be described by the random-field Ising model 
(RFIM) [2,3]: 



E{a m } 



CmO'm+l ) 



Ba„ 



(3) 



The Ising variable a m = ±1 corresponds to the two possi- 
ble values of the lattice dimerization Ai a t = ±Ao between 
neighboring kinks. Thus, the first term in Eq. (^) is the 
energy cost related to the occurence of kinks, where [i is 
the kink creation energy. As the RFIM Eq.(^) is an effec- 
tive model, obtained by integrating out small lattice fluctua- 
tions, jj, weakly depends on temperature and is renormalized 
by electron-electron interactions [2]. Comparing Eq. (|^) to 
Eq. (j2|) reveals that the role of the disorder is taken over 
by the random "magnetic" field h m , which we assume to 
have a Gaussian distribution with zero mean and correlator 
with the disorder strength e <C \i . 
Finally, the third term in Eq.(^|) describes the inter- 
chain interactions, which tend to establish a coherence be- 
tween the phases of the order parameter on different chains. 
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These interactions are taken into account in the chain mean- 
field approximation, so that the homogeneous "magnetic" 
field B is proportional to the average order parameter: 
B = W((a)). Here, the double brackets denote both the 
thermal and the random-field average. The validity of the 
chain mean-field approximation requires the interchain in- 
teraction energy to be sufficiently weak, W -C M, which 
holds in quasi-one-dimensional materials such as conjugated 
polymers. 
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Fig. 1. Phase diagram of the RFIM. 

We have analytically solved the model Eq. (^|) by notic- 
ing that the disorder averaged free energy of its contin- 
uum version has the form of a matrix element of the Green 
function describing the relaxation of a spin 1/2 in a time- 
dependent magnetic field. The coordinate along the chain 
plays the role of the (imaginary) time in which the relax- 
ation takes place. The Green function was found by solving 
the corresponding Fokker-Planck equation [2,3]. 

For a single chain (B tx W — 0) the ground state den- 
sity of solitons is then found to be proportional to the dis- 
order strength, 



e/n 



(4) 



The solitons occur randomly positioned along the chain and 
destroy the long-range order in the lattice dimerization at 
any temperature [2-4]. This turns out to be the same for 
interacting chains (W 7^ 0), as long as the disorder strength 
is larger than the critical value e c = 2W/J./3. However, 
for e < e c , the long-range bond order (LRBO) in the lat- 
tice dimerization is re-established as solitons are bound into 
pairs by the interchain interactions (soliton confinement). 
As a consequence, the density of neutral soliton pairs n v is 
found to be exponentially suppressed at low temperatures: 



W 2 ( 
— eXP l" 



(5) 



In Fig. 1 we plot the phase diagram for the existence of 
LRBO in the RFIM Eq. (^) using parameters typical for 
irans-polyacetylene. The result of a numerical calculation 
(stars) and our analytical result (solid line) are seen to be 



in excellent agreement. 

3. Magnetic Susceptibility in the Ordered Phase 

As neutral solitons carry a spin 1/2, they contribute 
to the polymer's magnetic properties. For e <C e c , solitons 
only occur in isolated pairs. Within the RFIM Eq. (|^), we 
have analytically calculated the pair size distribution p(R), 
which turns out to be sharply peaked around a typical pair 
size R = R* . As the antiferromagnetic exchange between 
the neutral solitons of a pair is a known (exponentially de- 
caying) function J(R) of the pair size [6], knowledge of p(R) 
allows us to calculate the magnetic susceptibility due to the 
disorder- induced solitons. 

At temperatures T larger than the typical exchange 
value J* = J(R*), we find the magnetic susceptibility 
to obey the Curie law, x(T) tx 1/T. At temperatures 
T <C J(R*), however, most of the spin pairs are in the 
singlet state, and the magnetic susceptibility (up to loga- 
rithmic corrections) reads [5]: 



X{T) 



(6) 



with a oc W 2 /e. Thus, the interplay between interchain 
interactions (W) and disorder (e) determines the low- 
temperature behaviour of the magnetic susceptibility. 

4. Conclusions 

Off-diagonal disorder induces neutral solitons in the 
ground state of trans-polyacetylene. These solitons con- 
tribute to the magnetic susceptibility. At low temperatures 
their magnetic response shows deviations from Curie behav- 
ior, and the thermal behavior of the susceptibility strongly 
depends on disorder and interchain interactions [Eq. (h)]. 
Our theory may explain the observed deviation from Curie 
behavior in Durham polyacetylene below 30K [7]. To get 
better insight into the parameter a, experiments should be 
extended to temperatures of the order of IK. 
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